Interacting random Dirac fermions in superconducting cuprates 
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We study the effects of quasiparticle interactions on disorder-induced localization of Dirac-like 
nodal excitations in superconducting high- Tc cuprates. As suggested by the experimental ARPES 
and terahertz conductivity data in Bi2Sr2CaCu20s+s , we focus on the interactions mediated by 
the order parameter fluctuations near an incipient second pairing transition d — > d + is. We find 
interaction corrections to the density of states, specific heat, and conductivity as well as phase 
and energy relaxation rates and assess the applicability of the recent localization scenarios for non- 
interacting random Dirac fermions to the cuprates. 



In the past few years, the problem of disordered two- 
dimensional (2D) Dirac fermions received much of atten- 
tion, as it provides an effective description for the ran- 
dom bond Ising model, network models of Quantum Hall 
plateau transitions, and some other statistical problems. 
Also, dynamical Dirac fermions can be used to conve- 
niently describe low-energy excitations in a variety of 
correlated systems, such as p-wave (e.g., Ru2Sr04) and 
layered d-wave (high- Tc cuprates) superconductors and 
superfluids {He^), and zero-gap semiconductors (e.g., 
graphene sheets). 

The Hamiltonian of a generic disordered gapless super- 
conductor and, in particular, planar d-wave system, pos- 
sesses an additional discrete symmetry of charge conjuga- 
tion which gives rise to as many as seven novel random 
Gaussian ensembles corresponding to different patterns 
of spin rotational and time reversal symmetry break- 
ing Q. Moreover, in a stark contrast with the conven- 
tional case of a normal metal with extended Fermi surface 
the density of states (DOS) of the non-interacting Dirac 
fermions in a 2D d-wave superconductor is strongly af- 
fected by disorder ||]. Furthermore, depending on the 
concrete model for disorder, such as isotropic versus pre- 
dominantly forward potential impurity scattering, the 
DOS of the random Dirac fermions can exhibit different 
low-energy asymptotic behaviors even within the same 
random ensemble 

The multitude of different regimes and crossovers pre- 
dicted for the non-interacting random Dirac fermions 
raises a question about their observability in such re- 
alistic d-wave systems as the high- Tc superconductors 
where quasiparticle interactions are believed to be im- 
portant. Moreover, the localization theory still remains 
incomplete without an extra input in the form of quasi- 
particle dephasing rate which controls the magnitude of 
disorder-induced localization corrections in the infinite 
system. 

Thus far, despite the continuing progress in under- 
staning of the non-interacting (de)localization phenom- 
ena in the d-wave systems, the above issues did not re- 
ceive enough attention. In the present paper, we fill in 



this gap by investigating the effects of physically relevant 
quasiparticle interactions on the localization properties of 
the Dirac-like nodal excitations in the superconducting 
cuprates. 

In the Nambu spinor representation, the quasiparticle 
(retarded) Green function reads as 



G^(e) = [(e - S^(6))ro + + Akfi]" 



(1) 



where we introduced the 2x2 unity matrix tq in addition 
to the triplet of the Pauli matrices fi_2,3. 

The bare quasiparticle spectrum which is composed 
of the normal state dispersion and the d^2_y2- 
symmetrical gap Ak becomes linear i?k = [^k + Ak]^^^ ~ 
[(v/(5k)^ -I- (vgSk)^]^/^ for momenta near each of the 
four gap nodes located at ±Ki^2 = ±(±fci?, A:f)/\/2. 
The pairs of orthogonal vectors v/ = (9^k/9k) and 
Vg — (3Ak/(9k) correspond to the components of the 
quasiparticle group velocity which are normal and tan- 
gent to the fudicial Fermi surface, respectively. 

By using the standard DOS definition 
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and linearizing the spectrum near the nodes, one readily 
recovers the linear DOS z/(e) = \e\/{TTv fVg) of the nodal 
low-energy excitations in a clean 2D d-wave supercon- 
ductor. 

The origin of the quasiparticle interactions in the 
cuprates remains a subject of an ongoing debate. The re- 
cent ARPES jl^ and optical (terahertz) conductivity |Tl| ] 
data in Bi2Sr2CaCu20^-\-s seem to rule out any short- 
ranged coupling that would have resulted in a tem- 
perature dependence of the inverse inelastic quasiparticle 
lifetime |12|. Instead, these data show an approximately 
linear behavior which is suggestive of the possibility of 
some sort of a quantum-critical behavior p^ ] . 

In the common scenario of a quantum-critical point 
(QCP), the long-ranged interactions between the quasi- 
particles result from their exchange by fluctuations of the 
new order parameter. The recent analysis of different in- 
cipient orderings singled out one of the second pairing 
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transitions dj.2_ „2 — » d^-2_y-2 +is(id'^y) as the most plau- 
sible candidate jl^ . Alternatively, one can also consider 
the long-ranged interaction due to charge fluctuations in 
the vicinity of a crystal structural instability Yet 
another type of interaction mediated by the exchange 
of antiferromagnetic fluctuations with momenta q «:! Q 
where Q = (±7r, tt) has been studied extensively in the 
context of the normal state of the cuprates and recently 
extended to the superconducting regime , albeit with 
no disorder included. 

In order to put all of the above quasiparticle interac- 
tions into one unifying framework we consider a generic 
form of the bosonic collective mode propagator 



Ti 



(l^/c)2 - (Q - q)2 - A2 - ViUfiLJ, q) 



(3) 



which describes the interactions in the spin, pairing, and 
charge channels for i = 0, 2, and 3, respectively. In Eq.(|^) 
the parameter controls the proximity to the QCP, and 
the fermion polarization functions n^(a;,q) account for 
the feedback effect of the nodal Dirac fermions on the 
spectrum of the collective mode. 

As regards modeling the disorder, we consider the 
experimentally relevant case of time-reversal and spin- 
rotational invariant scattering due to non-magnetic im- 
purities (class CI in the classification chartered in Ref. 
0) which we treat as scatterers with density rii cou- 
pled to the nodal Dirac fermions via the vertex C/k,k''r3- 
The latter accounts for arbirary strength and momen- 
tum dependence of disorder scattering and gives rise to 
the scattering T-matrix 



k'' 



C/k,k"f3G^»(e)]"'t/k",k'f3 (4) 



where the first and second terms in the brackets cor- 
respond to the Gaussian and Unitary limits, respec- 
tively 1^. In the earlier studies employing non-abelian 
bosonization and nonlinear a-model the lat- 

ter regime describing strong impurities has remained 
unattainable, and a markedly different behavior pre- 
dicted in this case [^,|| could not be readily addressed. 

The first insight into the random Dirac problem can be 
obtained from the self-consistent equation derived in the 
standard non-crossed diagrammatic approximation pi 
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rk,k'(e)G^Ke)7^k',k(' 



(5) 



As shown below, our neglecting the crossed diagrams is 
well justified by the large anisotropy of the Dirac spec- 
trum {vp/vg ^ 20) observed in the cuprates 

Provided that fii|Tk,k'(e)|/w/Wg << 1, for all the rele- 
vant energies and momenta, the solution of Eq. man- 
ifests the emergence of a new energy scale — /toE^(O) = 
7 << A which separates between the ballistic and diffu- 
sive regimes and gives rise to the finite DOS 
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n-^VfVg 7 



(6) 



In the non-crossed approximation, the spin and thermal 
conductivities which obey the Wicdcmann-Franz law ap- 
pear to be independent of disorder {h = ks = 1): 
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Att^T 47r2 



(7) 



while the charge conductivity receives non-universal ver- 
tex corrections |^. This is consistent with the fact that, 
unlike the charge of quasiparticles, their spin and energy 
are conserved, and, therefore, in the absence of localiza- 
tion they can both propagate diffusively Q . 

In the diffusive (e < 7) regime, the mean field DOS (|^) 
and the "universal" conductivities alike become sub- 
ject to further corrections which stem from both disorder- 
induced weak localization and interference between dis- 
order scattering and quasiparticle interactions. In order 
to compute these corrections one needs explicit expres- 
sions for all the gapless diffusion (V) and Cooperon (C) 
modes of the random Dirac problem beyond the ergodic 
limit studied in Ref. 

The abovementioned strong anisotropy of the Dirac 
spectrum in the cuprates implies as S> 1 which, apart 
from justifying the use of Eq.5 for calculating the 
fermion self-energy, allows one to resort to the stan- 
dard ladder approximation for the propagators 'D{C) = 
'^ijT^ij{Cij)Ti ® Tj expanded in the basis of the tensor 
products of the fi-matrices. 

By analogy with a normal metal, we first consider the 
ladders formed by one retarded and another advanced 
Green functions (RA-ladder) and find the following sin- 
gular contributions to the propagators of the soft modes 



I?„(Cy)(e,e',q) 
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TTVo Dq^ — i{e — e') 



(8) 



with the amplitudes rff = (1,1,1,1) and ?7p = 
(1,1,-1,1) corresponding to the T) and C propagators, 
respectively |l6|| . 

As previously pointed out the diffusion coefficient 
D appearing in Eq. (^) must satisfy the Einstein relation 
as — DiyQ/A. In our approach, this happens naturally, 
once the real part of the fermion self-energy ReS^(aj) 
from Eq.(|^) is accounted for in the ladder equations. 

As opposed to the case of a normal metal, in a super- 
conductor the gapless poles also appear in the RR (AA)- 
ladders due to the combination of processes of impurity 
scattering and Andreev refiection The correspond- 
ing propagators I?(C) are related to (||) by virtue of the 
charge conjugation performed on one of the two lines of 
the ladder: 



I?(C)(e,e',q) = - 5]C.,(2?y)(e, -e', q)(r2f*f2) < 



(9) 
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which can be cast in the form of Eq.(||) with ryf = 
(-1,1,-1,1) and 7?f = (-1,1,1,1). 

A straightforward analysis shows that it is the RR- 
Cooperon C which appears to be solely responsible for 
the first order weak localization DOS correction H 



1 



TT^D 
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(10) 



Eq.(|10D indicates that perturbation expansion remains 
well under control for energies above the characteristic 
scale ~ 7 exp(— 47r^(Ts) where it eventually breaks down. 

Besides the above soft modes, in the case of a nearly 
nested (square-like) normal state Fermi surface there 
might also exist additional (pseudo)Goldstone ones with 
gaps of the order of chemical potential fj, of the bare elec- 
trons. In the d-wave superconducting state, the corre- 
sponding poles occur in both RA- and RR(AA)-ladders 
when the sum of the two momenta carried by the lines 
equals Q, provided that the normal state dispersion sat- 
isfies « ~Ck±Q- In the Unitary limit, the presence 
of these extra modes reverses the sign of the localiza- 
tion correction ( |ic| ) in the energy interval < e < 7 
[^6| , resulting in a behavior which is reminiscent of that 
predicted in Refs. |^,^ for /i = 0. Nonetheless, at yet 
lower energies the DOS will eventually get suppressed, in 
accord with the scenario of Refs. [||j5[0]. 

Next, we calculate the interaction (Altshuler-Aronov 
type) corrections which stem from the interplay between 
disorder and quasiparticle interactions (||). As usual, the 
latter can be divided into the exchange and Hartree con- 
tributions. 

Upon computing the diffusion-dressed interaction ver- 
tices, we find that away from half-filling (2Ki_2 7^ Q) 
no diffusion poles occur for any finite transferred mo- 
menta Q 7^ which limits the subsequent discussion to 
the intra-node inelastic scattering. 

Furthermore, the vertices fi_3 undergo no singular dif- 
fusion dressing, as they correspond to the two spatial 
components of the Dirac fermion current operator. By 
contrast, both vertices fo and T2 which represent the 
Dirac fermion density and mass operators do develop 
diffusion poles in the RA- and RR(AA)-channels, respec- 
tively. 

In light of the above, from now on we focus on the 
coupling V^'{u!, q) mediated by the fluctuations of a sec- 
ondary pairing order parameter with Q = which was 
suggested as a possible scattering mechanism in the high 
temperature (ballistic) regime ]T^ . Here we restrict our 
attention to the case of the is secondary pairing, for the 
case of the idxy pairing appears to be somewhat more 
intricate ]T7[ |. Unlike any other before mentioned types 
of interactions, it is strongly enhanced by disorder, thus 
producing the diffusive DOS correction 
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In the conserving approximation, the above vertex renor- 
malization is necessarily accompanied by the singular po- 
larization of the disordered Dirac fermions 



n^(cj,q) = I'D In 
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Dq^ 



(12) 



Close to the QCP (A^ < V2V0), the term (|T|) domi- 
nates over the other terms in the denominator of Eq. (||) , 
thereby resulting in the intrinsically attractive coupling 
V^{uj,q) « — l/n2'(w,q) which yields a positive DOS 
correction 



^^Inlln^l 



(13) 



Thus, in contrast to the exchange DOS correction in a 
normal metal, in our case 5vex appears to be substan- 
tially weaker than the effect of weak localization ([l0|). 

Farther away from the QCP (A^ > V2Vq), the bosonic 
propagator (|^) becomes effectively short-ranged, and the 
(positive) interaction correction to DOS becomes explic- 
itly dependent upon the coupling strength: 5f,xv{e) ~ 
(VJji'o/A^) ln7/e, while the corresponding Hartree term 
vanishes identically. 

The latter regime appears to be partly similar to the 
case of a short-ranged non-singular ferromagnetic cou- 
pling of the local quasiparticle spin densities in a gapless 
bulk superconductor which yields the interaction DOS 
correction of the same functional form and sign as the 
localization one ||l8| . 

The DOS correction (|l|) can also be deduced from the 
singular contributions to thermodynamic quantities, e.g., 
specific heat 



tan ^ 



/mHf^(w, q) 
iMI?(^ 
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^^Inlln^l 



(14) 



calculated for A 0. Accordingly, the exchange con- 
tribution to the conductivity correction is given by the 
expression 



(5ea;Crs(T) 
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(15) 
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which also shows a much weaker temperature depen- 
dence than the logarithmic weak-localization correction 
5wi(js{T) — — ln(7T0(T))/27r^ controlled by the inelastic 
phase-breaking time (T) Q . 

The latter can be estimated as a time interval for which 
the accumulated phase uncertainty 



q 



did , / N 



1 — exp(icji — Dq^t) 



(16) 



becomes of order unity [gOj . From the condition 5'^{t^) 
1, we obtain the estimate 



T 



IGttcts In 7/r 



(17) 



that should be contrasted against the apparent phase- 
breaking rate oc T^/^ found to describe both, experimen- 
tally and theoretically, quasiparticle transport in the nor- 
mal state of the cuprates pT[|. 

Also, the dephasing time(0) differs from the inelastic 
energy relaxation (or, "out-scattering") time [19| 



duj , , , w s , , a; , , 
-[coth(-)-tanh(-)] 



5]/my2^(w,q) 



1 



(18) 



main predictions of the non-interacting localization sce- 
nario of Refs. Lastly, we provided the so far miss- 
ing ingredient of the localization theory by computing 
the inelastic dephasing and energy relaxation. 

We note, in passing, that the latter does not ac- 
count for the possibility of the formation of a strongly 
anisotropic network of delocalized resonant impurity 
states, as proposed in Ref. |2^. However, as far as the 
thermodynamic and transport properties of the rest of 
the electronic spectrum are concerned, our findings pro- 
vide further support for the attempts to find manifesta- 
tions of the quasiparticle localization-related phenomena 
in the superconducting cuprates |^3|] . 

Meanwhile, the high precision heat transport measure- 
ments [p4| , p5| might have already provided such an impor- 
tant evidence of weak localization as positive magneto- 
(thermal)conductivity. Should this interpretation of the 
data of Ref. proposed in prove to be correct, the 
phase relaxation time ( [T7| ) may as well turn out to be 
experimentally accessible. 

The authors are grateful to A. Altland for his valuable 
comments on the case of the id secondary pairing and 
also on Eq.(ll). This research was supported by the NSF 
under Grants DMR-0071362 (DVK) and DMR-9703388 
(AGY), by RFBR (IVG), and partly by INTAS. 



For T << 7/0-5 the solution of this self-consistent equa- 
tion can be found in the form 



167r CTg In 7/r 



(19) 



Unlike which is determined by the processes with en- 
ergy transters ui ^ '''<^^' ''"e receives contributions from 
the frequencies up to w ~ T. In the ballistic (T > 7) 
regime, t~^(T) becomes a linear function of tempera- 
ture, in agreement with the earlier theoretical ||l3|| and 
experimental (l^Jll[l results. 

In conclusion, we carried out the analysis of the effects 
of the experimentally relevant quasiparticle interactions 
mediated by fluctuations of a secondary pairing order pa- 
rameter (is) on the localization properties of the nodal 
quasiparticles in disordered high- Tc superconducors. In 
the course of this study we revisited the problem of local- 
ization of noninteracting random Dirac fermions and ex- 
plicitly solved for all the relevant diffusion and Cooperon 
modes in a 2D d-wave superconductor with arbitrarily 
strong impurities. We found that, despite being singu- 
lar as a function of transferred energy and momentum, 
the interaction in question generates only subdominant, 
albeit positive, DOS, specific heat, and conductivity cor- 
rections and, therefore, does not necessarily alter the 
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